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The present B —> nK data is studied in the context of the standard model (SM) and with new 
physics (NP). We confirm that the SM has difficulties explaining the B — > nK measurements. By 
adopting an effective-lagrangian parametrization of NP effects, we are able to rule out several classes 
of NP. Our model-independent analysis shows that the B ttK data can be accommodated by NP 
in the electroweak penguin sector. 



PACS numbers: 13.25.Hw, ll.30.Er, 12.15Ff, 14.40.Nd 

The i3-factories BaBar and Belle have measured (most 
of) the branching ratios and CP asymmetries for the var- 
ious B — > TTTT and B — > ttK decays, and these can be used 
to search for physics beyond the Standard Model (SM). 
By using flavor SU(3) symmetry to relate these processes 
[IIIELS, 4j_^,Je,7i, several analyses were able to constrain 
the SM parameters, and to look for signs of New Physics 
(NP). The advantage of this approach is that one takes 
into account a large number of processes. The disadvan- 
tage is that one has to deal with unknown effects related 
to the breaking of SU(3) symmetry. Also, B ^ tttt de- 
cays involve the quark- level processes b — > dqq (q = u,d), 
while B ttK receives contributions from b sqq. If 
there is NP, it could affect b ^ d and b ^ s processes 
differently. 

For this reason, there are advantages to considering 
B ttK decays alone. As we will see, these pro- 
cesses contain enough information to constrain the SM 
parameters. Within the diagrammatic approach the 
amplitudes for the four B ttK decays can be writ- 
ten in terms of seven diagrams: the color-favored and 
color-suppressed tree amplitudes T' and C", the gluonic 
penguin amplitudes P' and P^^, the color- favored and 
color-suppressed electroweak penguin amplitudes P^„,- 
and -Pew-i the annihilation amplitude A'. (The 

primes on the amplitudes indicate b —f s transitions.) 

In Ref. [3, the relative sizes of the amplitudes were 
roughly estimated as 
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where A ^ 0.2. These estimates are expected to hold 
approximately in the SM. Thus, to 0(A), we can ignore 



all diagrams but P', T' and P^„, in our B ttK ampli- 
tudes. We will perform a fit of the present B ttK data 
- the goodness or badness of the fit should not be much 
affected by the inclusion of the smaller amplitudes. 
The four amplitudes can then be written as 
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where we have explicitly written the dependence on the 
weak phase (including the minus sign from V^^Vt^ [P']), 
while the amplitudes contain strong phases. (The phase 
information in the Cabibbo-Kobayashi-Maskawa (CKM) 
quark mixing matrix is conventionally parametrized in 
terms of the unitarity triangle, in which the interior (CP- 
violating) angles are known as a, /3 and 7 0-) 

We have one additional piece of information: within 
the SM, to a good approximation, the diagram P^,^ can 
be related to T' using flavor SU(3) [l^: 
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Here, the Ci are Wilson coefficients ^^'^ 
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With the above relation, the B ttK observables 
contain five theoretical parameters: |P'|, |T'|, /3, 7, and 
one relative strong phase, 5. The phase (3 can be taken 
from the measurement of sin 2/3 in B^{t) — > J/ipKs'. 
sin 2/3 = 0.726 ± 0.037 |T3, leaving four theoretical un- 
knowns. However, there are a total of nine B — > ttK mea- 
surements: four CP-averaged branching ratios and five 
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TABLE I: Branching ratios, direct CP asymmetries A^ir, and 
mixing-induced CP asymmetry Aindir (if applicable) for the 
four B nK decay modes. 



A 



-+ 



-T'p^^ ^ P' - P' p*^ P"^ 





Mode 


Bi?(10-^) 




B+ 




24.1 ± 1.3 


-0.020 ± 0.034 


B+ 




12.1 ± 0.8 


0.04 ±0.04 


B'l 


■K-K+ 


18.2 ± 0.8 


-0.109 ±0.019 






11.5 ± 1.0 


-0.09 ±0.14 0.34 ±0.28 



CP asymmetries (Table |T| [13). Within the parametriza- 
tion of Eq. three of these are independent of the 
theoretical parameters: the direct CP asymmetries in 
B+ TT+K^ and B° Tr°K° are predicted to vanish, 
and the indirect CP asymmetry in B'^ — > tt'^K'^ measures 
sin 2(3. The remaining six observables are functions of 
the four theoretical parameters, so we can perform a fit 
to obtain these quantities. 

Using the parametrization of Eq. (O for the B — > ttK 
amphtudes, we find that x^^^/d.o.J. = 15.6/5 (0.8%), 
indicating an extremely poor fit. (The number in paren- 
theses indicates the quality of the fit, and depends on 
Xmin ^i^d d.o.f. individually. 50% is a good fit; fits which 
are substantially less than 50% are poor.) This is not a 
new result - other analyses have made a similar obser- 
vation P, IE S S IE 0- It shows that present data is 
inconsistent with the naive implementation of the SM. 
Our parametrization is therefore incomplete. There are 
two ways to make modifications. Either we work within 
the SM, or we add new physics. We address these possi- 
bilities in tm'n. 

We begin with the SM, but abandon the relation be- 
tween P'J^^y and T' [Eq. Q]. We now have six theoretical 
parameters: \T'\, |P^„-|, 7, and two relative strong 

phases. In this case, the fit is good: x^^^/d.o.f. — 2.7/3 
(44%). The fit also gives a central value of 7 = 59°. This 
is consistent with the value of 7 obtained via a fit to 
independent measurements: 7 — 62^}2 01- (Because 
these errors are not gaussian, we do not include this in- 
formation in our fit at this stage.) On the other hand, 
the fit also gives \Pl,,jT'\ = 1.55, which is far from its 
SM value of 0.65 ± 0.15 0|. Thus, while it is possible 
to explain the present B -kK data by treating P'j^^- 
and T' independently, it is difficult to understand how 
\P'e„/T'\ could be 6cr away from its SM value. 

The second modification is to take into account the 
smaller (neglected) amplitudes. Including the O(A^) di- 
agrams, the B — > ttK amplitudes take the form 
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In this case, P^^y is not independent of the amplitudes 
T' and C We have 

P^^^. = 3 cg + cio^^y, ^ C')±?^^^^i?(r' - C) , 

4 Cl + C2 4 Cl — C2 

p-^ ^ 3 C9±CK,^ ^ ^,y^C_^-C^ 
4 Cl + C2 4 Cl - C2 

With these relations, we now have eight theoretical 
parameters: \P'\, IP^^Ii l^'l' 7, and three relative 
strong phases. With nine pieces of experimental data, we 
can still perform a fit, which is acceptable: Xmm/'^-'^-f- = 
0.7/1 (40%). In addition, we find 7 = 64°, consistent 
with independent measurements. However, the fit gives 
\C' /T'\ — 1.8, which is far larger than naive estimates. 
Other analyses have also found the C must be very big 
to explain the B ^ ttK data H H |1 0- In Ref. 0, it 
is argued that final-state interactions (FSI) can increase 
the size of C". However, in that case, the authors were 
attempting to explain \C' /T'\ ~ 0.5 (which comes from 
the joint fit to £? ^ ttK and B — > tttt decays |3|)- Even 
with FSI, it is difficult to see how C can be increased to 
about twice as large as T'. 

It is therefore extremely difficult to explain the current 
B ttK data within the SM alone. Instead, one must 
consider the effect of new-physics operators. These can 
be included in i3-physics analyses in a model-independent 
way . We suppose that there are NP contributions to 
b sqq transitions which are roughly the same size as 
the SM 6 — > s penguin operators. The NP contributions 
take the form Onp ~ sTibqTjq (g = u, d, s, c), where the 
Tij represent Lorentz structures, and colour indices are 
suppressed. There are a total of 20 possible NP opera- 
tors, each of which can in principle have a different weak 
phase. The NP contributes to the decay B ^ f through 
its matrix elements (/| O^p \B), which can be written as 



{f\0''^',^\B)^Ake''f'le''l 



(7) 
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where 0| and 5^ are the NP weak and strong phases 
associated with the individual matrix elements. How- 
ever, the key point is that the NP strong phases are 
very small. The reasoning goes as follows. Strong phases 
arise from rescattering. In the SM, the (large) tree dia- 
gram (T') b — > sec can rescatter into the c-quark penguin 
P^, possibly giving it a strong phase of 0(1). Note that 
\Pc/T"\ < 10%. That is, in the SM the diagram respon- 
sible for the rescattering is considerably larger than the 
diagram which receives the strong phase. On the other 
hand, the NP rescattering can only come from the NP 
matrix elements themselves. Assuming the same sup- 
pression factor, the NP strong phases are 0(10%), which 
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is negligible, to a good approximation. Note that this is 
a quite general result and applies to all NP models. 

The neglect of NP strong phases allows for a great sim- 
plification. For a given type of transition, all NP matrix 
elements can now be combined into a single NP ampli- 
tude, with a single weak phase: 

(/I O^'^! \B) = ^^e^*' , (8) 

where q — u,d, s, c. (Throughout the paper, the symbols 
A and $ denote the NP amplitudes and weak phases, 
respectively.) B — > ttK decays involve only NP param- 
eters related to the quarks u and d. These operators 
come in two classes, differing in their colour structure: 
Sa^ibaQf^Tjqf:! and SaTibpqpTjQa {q = u,d). The ma- 
trix elements of these operators can be combined into 
single NP anmlitudes, denoted A'^'^e'^'" and y^'^-'e'*^"', 
respectively [l6j- Here, and are the NP weak 
phases; the strong phases are zero. Each of these con- 
tributes differently to the various B ttK decays. In 
general, ^ and ^ <^'^ . Note that, despite 

the "colour-suppressed" index C, the matrix elements 
^/G,9gi*g g^j,g j^Q^ necessarily smaller than the ^''^e**? . 

The B — i- ttK amplitudes can now be written in terms 
of the SM amphtudes to 0(A) [P^w ^^'^ related 
as in Eq. 0], along with the NP matrix elements [l^ : 

A+o = -P' + A'^-'^e'^'^^ , 
y^A^+ = P' - r e'^ ^ P'^,,, 

_j_ j^l ,comb ^i^' a! e.^^'^ 

A~+ =P' -T' e'^ - ^"=^'"6**"" , 
V2A°° = -P' - P^„, + ^'.™"''e'*' + ^'^'''e'*^^ , (9) 

where ^'-co^bg**' ee ~y^''«e**'" + A' e^'^'-^ . There are now 
a total of 11 theoretical parameters: \P'\, \T'\, |^''™™''|, 
\A'°-''^l \A'^''^\, 7, 3 NP weak phases and two relative 
strong phases. With only 9 experimental measurements, 
it is not possible to perform a fit. It is necessary to make 
some theoretical assumptions. 

We assume that a single NP amplitude dominates. 
There are an infinite number of choices, but we consider 
the following four possibilities: (i) only A'''^°"^^ ^ 0, (ii) 
only ^"^'^ 7^ 0, (ui) only A"'''^ ^ 0, (iv) yl'^-^e^*-'' = 
_4/c,dgi<[.f ^ ji,comb ^ Q (isospin-conserving NP). 

In the first three cases there are seven parameters: 
three amplitude magnitudes, 7, one weak NP phase and 
two relative strong phases. However, for the type of NP 
characterizing the fourth fit, all B ttK amplitudes and 
their CP-conjugates contain two combinations of ampli- 
tudes. These can be written as follows: 

p^^giSjvPgi'I'lvp ^ —P' + A"^'''^e^'^''^ 

P^pe^'^^'-e-**"^ = -P' + ^"^•''e-'*^'' , (10) 

with Pjvp ^ Pjvp- However, note that the real parts 
of these quantities are equal: P^p cos((5„p -t- ^np) = 



Pjvp cos(5jvp — $ivp)- Thus, one variable, say S^p, can 
be written as a function of the other three. That is, in 
this case there is one fewer degree of freedom, and the 
B ttK amplitudes contain six unknown parameters: 
Pjvp, Pjvp, T', 7i ^np and St'- 

In cases (i) , (iii) and (iv) , there are more measurements 
than unknowns, and we can perform a fit. On the other 
hand, parametrization (ii) makes the same predictions as 
the SM: Ad^r{B+ Tr+ K°) ^ AdUB^ Tr°K°) = 
and Aindir{B'^ tt'^K'^) = sin 2/3. Thus, in this case, 
there are only six observables, and we cannot determine 
the seven theoretical parameters. 

Using Table our results are: (i) x^j„/d.o./. = 
1.9/2 (39%), (iii) xLJd.o.f. = 9 A/2 (0.9%), (iv) 
xl^iJd.o.f. = 3.9/3 (27%). Thus, based on the fit qual- 
ity only, we conclude that fit (i) is acceptable, fit (iv) is 
marginal, and fit (iii) is poor. 

However, these fits also give values for the CP angle 
7: (i) 7 = 64.2°, (in) 7 = 31.8°, (iv) 7 = 37.8°. These 
can be compared with the value obtained from a fit to 
independent data, 7 = 62^12° Note that this latter 
value of 7 includes limits on B'^-B^ mixing. However, 
the NP considered here will, in general, also lead to ef- 
fects in this mixing. Thus, technically, in considering this 
type of NP, the B^-B'^ mixing data should be removed 
from the 7 fit. In practice, though, this will not make 
much difference. We therefore continue to use the best- 
fit values of 7 with 7 = 62^i[2 as the independent value, 
but the reader should keep this caveat in mind. 

Note also that any explanation of the B ttK data 
using new physics must also reproduce the SM value of 7. 
This demonstrates that, in looking for NP, it is important 
to use all handles available, and not simply concentrate 
on measurements of the CP phases. 

We incorporate the information on 7 by adding a con- 
straint to the data, so that we now fit to all the B ttK 
data and 7 = 62 ±11°. (Note that the 7 constraint 
is not a true experimental number - it has some the- 
oretical input - and so its inclusion in the fit must be 
viewed with some prudence.) With this added input, we 
can now perform a fit in parametrization (ii). We find: 

(i) xLJd-o.f. = 1.9/3 (59%), (ii) xLJd-o.f. - 2.7/3 
(44%), (iii) xLJd.o.f. = 9.4/3 (2%), (iv) xLJd.o.f. = 
6.7/4 (15%). We conclude that fits (i) and (ii) are good, 
while fit (iv) is poor, and fit (iii) is ruled out. We do not 
consider fit (iii) further. 

However, we have still not included all the information 
at our disposal. In the fits, we find that (i) St' = —58.4°, 

(ii) St' = -26.2° or 68.8°, and (iv) St' = -47°. On 
the other hand, the diagram T' is governed by the CKM 
matrix elements V*i^Vus, and so its strong phase can arise 
only from self-rescattering. Thus, like the new-physics 
amplitudes, the strong phase of T', St', is expected to 
be very small. This requirement gives us an additional 
handle. We incorporate this by adding a constraint to 
the data: we require 6t' = 0± 10°. Since this is purely 
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theoretical, it is obviously not on the same footing as 
the B ttK data. However, since we only want to see 
if particular NP models give a good fit, it is sensible to 
include this information among the constraints. 

Including the constraint on St' , we find (i) 
xLJd-o.f. - 2.2/4 (70%), (ii) xLJd-o.f. = 5.9/4 
(21%), and (iv) xlnn/d-o.f. = 14.3/5 (1%). We therefore 
find that (i) is a good fit, (ii) is disfavored, and (iv) is a 
very poor fit. 

Of the four new-physics models examined in this pa- 
per, only one produces a good fit to the B ttK data 
and the various imposed constraints on 7 and St' ■ It 
is case (i), A'^"""^^ ^ 0. In this model, the best-fit val- 
ues of the theoretical parameters are \T' /P'\ = 0.22 (in 
line with theoretical expectations), |^'.™mfc/p/| ^ q gg^ 
= 100°, 5p' = -10°. We therefore find that the NP 
amplitude must be sizeable, with a large weak phase. 

This class of NP models essentially corresponds to a 
modification of the SM electroweak penguin amplitude, 
as explored in Refs. (H, ||3(, il^j . In Ref. (ll the weak phase 
of the electroweak penguin was modified, meaning that 
the NP operator is of the form {V — A) x {V — A). Here, 
we allow any form for the operator, so that this is a more 
general solution. NP models which can lead to this in- 
clude Z- and Z'-mediated flavour-changing neutral cur- 
rents |l8l Il9l| or supersymmetry with i?-parity breaking. 

Fit (ii) {A"^''^ ^ 0) is marginal, but not completely 
ruled out (it does give a value for \T' /P'\ which is about 
three times larger than expectations) . This is a NP solu- 
tion which has not been considered before. It can arise, 
for example, in supersymmetric models with i?-parity 
breaking. Fit (iii) {A"^'"^ ^ 0) yields a poor fit, and this 
class of NP models can therefore be ruled out. We also 
rule out isospin-conserving models of NP [fit (iv)]. These 
include new physics whose principal effect is to generate 
an anomalous gluonic quadrupole moment . 

A word of caution: one has to be careful about ruling 
out particular models of new physics. Any specific NP 
model will, in general, lead to more than one NP oper- 
ator, and the more general case can be used to explain 
the B — > ttK data. 

To summarize, we have presented a study of the cur- 
rent B ttK data. The standard model (SM) has 
great difficulty accounting for these measurements. De- 
pending on the parametrization, one obtains a poor fit, 
or values for the SM parameters which are greatly at 
odds with our present understanding. For models of new 
physics (NP), we adopt a model-independent, effective- 
lagrangian parametrization of the NP effects. There are 
three possible (complex) NP parameters which can affect 
B ^ ttK decays, denoted y^''™™*-, y^'c," and A"'^'^. We 
consider four classes of NP models: (i) only _4'>'=°™'' ^ Q, 
(ii) only ^ 0, (iii) only A'"^'' ^ 0, (iv) isospin- 

conserving NP: A'^-^e"^'^ = A'^-'^e''^'"'' , A' """^'' = 0. Of 
these, the classes of models (ii), (iii) and (iv) are ruled out 
or disfavored. Only model (i) explains the data satisfac- 



torily. It corresponds to a modification of the electroweak 
penguin (EWP) amplitude. Note that, while other stud- 
ies also consider specific models of NP in the EWP, our 
analysis is completely model independent. 
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